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Me" study  the  flow  of  two  immiscible  fluids  of  different  viscosities  and 
equal  density  through  a  pipe  under  a  pressure  gradient.  This  problem  has  a 
continuum  of  solutions  corresponding  to  arbitrarily  prescribed  interface 
shapes.  The  question  therefore  arises,  which  of  these  solutions  are  stable 

and  thus  observable.  Experiments  have  shown  a  tendency  for  the  thinner  fluid 

>  - 

to  encapsulate  the  thicker  one.  This  has  been  "" explained"  by  the  viscous 
dissipation  principle,  which  postulates  that  the  amount  of  viscous  dissipation 
is  minimized  for  a  given  flow  rate.  For  a  circular  pipe,  this  predicts  a 
concentric  configuration  with  the  more  viscous  fluid  located  at  the  core.  A 
linear  stability  analysis,  which  is  carried  out  numerically,  shows  that  while 
this  configuration  is  stable  when  the  more  viscous  fluid  occupies  most  of  the 
pipe,  it  is  not  stable  when  there  is  more  of  the  thin  fluid.  Therefore  the 
dissipation  principle  does  not  always  hold,  and  the  volume  ratio  is  a  crucial 
factor.  <* - - ' 
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SIGNIFICANCE  AND  EXPLANATION 


When  two  immiscible  fluids  of  different  viscosities  and  nearly  equal 
densities  flow  through  a  pipe  under  a  pressure  gradient,  experiments  have 
suggested  that  no  matter  what  the  initial  configuration,  the  less  viscous 
fluid  eventually  encapsulates  the  more  viscous  fluid.  This  property  has  been 
observed  for  both  low  and  high  Reynolds  number.  For  example,  in  the  pipeline 
transport  of  viscous  oils,  the  pressure  gradient  can  be  reduced  by  adding 
water  because  it  tends  to  coat  the  pipe  wall.  Another  example  arises  in  the 
spinning  of  bicomponent  fibers  such  as  nylons,  when  two  polymer  omits  are 
extruded  through  a  tube. 

Mathematically,  anti-planar  shear  flow  (exclusively  axial  flow  with  only 
one  non-zero  component  of  velocity  which  depends  on  the  coordinates 
perpendicular  to  the  axial  coordinate)  at  low  Reynolds  number  in  a  cylindrical 
pipe  of  arbitrary  cross-section  has  a  continuum  of  solutions:  for  any  pre¬ 
assigned  interface  shape  there  is  a  possible  flow.  Previous  theoretical 
attempts  to  account  for  the  unique  observed  flow  have  relied  on  the  'viscous 
dissipation  principle*  which  states  that  the  most  favored  configuration  is  the 
one  which  minimizes  viscous  energy  dissipation  for  given  volume  flux.  This 
holds  for  low  Reynolds  number  one-component  flow  but  has  no  mathematical  basis 
for  bicomponent  flows.  Rather,  we  thought  that  a  stability  analysis  ought 
to  be  done.  In  fact,  for  axisymmetric  flow  in  a  pipe  of  circular  cross- 
section,  linear  stability  results  for  long  waves  have  been  calculated 
numerically  by  Hickox  (1971).  However,  he  only  computed  for  the  case  in  which 
the  less  viscous  fluid  is  encapsulated  by  the  higher  viscosity  fluid.  He 
found  this  case  to  be  unstable  at  all  Reynolds  numbers.  He  did  not  consider 
the  case  In  which  the  more  viscous  fluid  is  centrally  located  to  test  for  the 
stability  predicted  by  the  viscous  dissipation  principle.  We  have  therefore 
reconsidered  the  linear  stability  problem  but  without  approximations.  The 
results  shows  that  the  viscous  dissipation  principle  does  not  hold  and 
stability  depends  primarily  on  the  ratio  of  the  radii  of  the  core  fluid  and 
the  pipe.  Some  of  our  results  show  qualitative  similarities  with  results  of 
Yih  (1967)  on  the  stability  of  plane  Couette  flow.  The  instability  is  due  to 
a  mode  which  is  neutrally  stable  when  the  viscosities  of  the  two  fluids  are 
equal.  It  is  also  neutrally  stable  in  the  limit  as  the  viscosities  became 
large  so  that  the  instability  is  purely  a  finite  Reynolds  number  effect.  This 
can  also  be  seen  in  Hickox 's  work  and  provides  a  reason  as  to  why  the  viscous 
dissipation  principle  is  not  always  in  agreement  with  analysis  of  stability. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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A  aajor  problem  in  the  theory  of  bicomponent  flows  lies  in  their  nonuniquenees:  The 
position  of  the  Interface  is  one  of  the  unknowns,  but  the  equations  of  motion  may  permit  an 
infinite  number  of  different  interface  configurations*  The  question  thus  arisess  which 
interface  positions  are  stable  and  thus  observable?  A  number  of  experiments  (1-3),  (51, 
(?),  (11-14]  have  shown  a  tendency  for  the  thinner  of  the  two  fluids  to  move  into  regions 
of  high  shear.  This  has  led  to  the  conjecture  that  the  stable  flow  can  be  characterized  by 
a  minimisation  principle  (2],  (61.  The  amount  of  viscous  energy  dissipation  should  be 
minimized  under  appropriate  constraints. 

In  this  paper,  we  focus  on  antiplane  shear  flow  in  a  cylindrical  pipe.  This  situation 
has  applications  in  the  transport  of  oil,  where  the  flow  rates  can  be  increased  by  adding 
som  (less  viscous)  water.  Another  application  arises  in  the  co-extrusion  of  two  polymers, 
e.g.  in  the  manufacture  of  nylon  fibers  (2).  The  Implications  of  the  "viscous  dissipation 
principle"  are  discussed  in  chapter  2.  The  flow  of  a  single  fluid  is  described  by  the 
variational  formulation 
(1.1) 


min  /  ^  u(Vu)2  -  Gu  . 

udnj(fl) 
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8m,  the  velocity  has  the  fora  u,  -  u(x,y)  •  -G  la  a  constant  pressure  gradient,  and 
9  Is  ths  cross-section  of  the  pips,  for  a  two-component  flow,  y  is  a  discontinuous 
function  of  x  and  y.  For  each  choice  of  y,  ths  functional  (1.1)  has  a  minimum, 
leading  to  a  solution  of  ths  Mavier-Stokes  equations.  Ms  lspose  ths  constraint  that  the 
volume  occupied  hy  each  fluid  is  given,  thus  y  is  allowed  to  vary  in  a  set  of  functions 
taking  constant  values  in  two  regions,  whose  measures  are  prescribed.  The  conjecture 
expressed  in  (2),  [6)  can  then  be  stated  as  minimising  (1.1)  not  only  with  respect  to  u, 
but  also  with  respect  to  y .  This  is  equivalent  to  saying  that  the  amount  of  viscous 
dissipation  (or  work,  since  the  two  are  equal)  is  minimised  for  given  flow  rate,  and 
maximised  for  given  pressure  gradient.  Alternatively,  we  may  say  that  the  pressure 
gradient  la  minimised  for  given  flow  rate,  or  the  flow  rate  maximised  for  given  preesure 
gradient. 

for  a  general  cross-section  of  the  pipe,  we  do  not  know  whether  this  variational 
problem  always  has  a  solution.  Me  show  that  the  existence  of  the  alnimlser  la  implied  by 
an  a  priori  estimate  on  the  length  of  the  interface  curve,  but  we  do  not  know  how  to  obtain 
such  an  estimate.  For  a  circular  pipe,  the  symmetry  leads  to  an  enormous  simplification. 
The  minimi ser  can  be  constructed  explicitly  by  elementary  means.  Xt  is  a  concentric 
configuration  with  the  sure  viscous  fluid  located  at  the  core. 

For  this  case  of  a  circular  pipe,  we  assess  the  validity  of  the  disslpetlon  principle 
by  performing  a  linear  stability  analysis.  This  is  done  in  chapter  3.  The  linearised 
equations  are  solved  maserlcally  by  Orssag’s  (8]  method  of  POurier-Chebyshev  expansions. 

The  linear  stability  of  bicomponent  flows  was  studisd  previously  by  Tlh  [14]  and  Hlckox  (4) 
using  perturbation  expansions  for  long  waves.  Tih  studies  plane  Oouette  and  plane 
Poiseuille  flow  with  a  single  flat  interface.  They  did  not  pose  the  problem  of 
selection.  In  these  cases  the  dissipation  principle  would  predict  instability  for 
Poiseuille  flow,  while  It  gives  no  prediction  for  Couette  flow.  Tlh  found  instability  in 
the  case  of  poiseuille  flow;  for  Oouette  flow,  the  stability  depends  on  the  volume  ratio  of 
the  two  fluids.  Hlckox  (4]  studied  concentric  flow  in  pipes  with  the  less  viscous  fluid 
located  at  the  core  and  finds  instability.  Re  does  not  consider  the  reverse  case  where  the 
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•on  vlacoua  fluid  occupies  the  cor*  as  prsdictsd  by  ths  dissipation  principle.  Both  Yih 
and  Hickox  find  that  th*  critical  eigenvalue  becomes  imaginary  for  s*ro  Rsynolds  niabtr, 
hsnes  th*  stability  or  instability  of  th*  interface  is  strictly  •  finite  Rsynolds  number 
effect,  this  shoes  that  encapsulation  Is  a  nonlinear  phenoamnon,  governed  by  th*  Havler- 
Stokes  equations  rather  than  Stokes  equations.  It  nay  be  noted  that  th*  Wavier- Stoke* 
equations  do  no  arise  froa  a  variational  principle. 

In  our  computations ,  ve  studied  stability  tor  both  cases  in  which  th*  less  and  nor* 
viscous  fluids  are  centrally  located.  The  forner  Is  always  unstable.  For  the  latter,  th* 
aost  Important  factor  was  found  to  be  the  volume  ratio  of  the  two  fluids.  When  th*  thicker 
fluid  occupies  aost  of  th*  pipe,  we  find  stability  as  the  viscous  dissipation  principle 
predicts.  However,  as  th*  radius  of  the  cor*  is  reduced  past  about  0 . 7r2 ,  where  R2  is 
th*  pip*  radius,  we  find  instability.  This  shows  that  in  any  event  th*  dissipation 
principle  is  not  strict.  It  is  perhaps  of  interest  that  th*  solution  of  design  problea 
(see  §10  of  [5] )  for  the  thickness  of  the  water  layer  which  will  naxials*  the  flux  of  oil 
in  th*  core  of  th*  pip*  le  in  th*  region  of  stability  for  sufficiently  long  waves.  We  find 
that  th*  critical  radius  ratio  increases  with  th*  wav*  nuaber  of  th*  perturbation.  In 
reality,  instability  of  short  waves  aay  be  suppressed  by  surface  tension,  which  was  not 
included  in  our  considerations. 

Th#  reader  aay  be  inclined  to  think  that  our  results  con fir*  the  dissipation  principle 
when  th*  voliaw  ratio  is  biased  towards  th*  aore  viscous  fluid,  we  feel,  however,  that 
great  caution  is  advisable  her*.  The  situation  we  have  dealt  with  has  a  rotational 
syaaetry,  and  aynswtric  solutions  to  syaastrlc  problea*  always  have  a  special  status.  It 
is  thus  quit*  natural  that  th*  concentric  configuration  should  be  preferred  over  others  in 
certain  situations,  whatever  th*  mechanism  aay  be. 

On  th*  otlw  ;  hand,  we  find  that  the  growth  rates  of  th*  unstable  aodea  for  th*  ainiaua 
dissipation  case  are  smaller  by  on*  or  two  orders  of  magnitude  than  for  th*  revere*  case  in 
which  th*  thin  fluid  is  at  th*  cor*.  Thus  th*  configuration  selected  by  th*  dissipation 
principle  appears  to  be  at  least  less  unstable  than  some  other  interface  positions. 
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2.  The  dissipation  g|aclgl| 

we  consider  flow  in  •  horizontal  cylindrical  pipe,  whoso  cross-soction  Q  is  s 
bounded  domain  with  s  (Booth  boundary.  The  pips  is  occupied  by  two  fluids  of  equal 
densities,  but  different  viscosities  u 1  and  u2  >  u^.  We  study  statlonsry  antiplane 
shear  flow,  i.e.,  if  x  and  y  denote  transverse,  and  s  the  longitudinal  coordinate, 
then  the  velocity  field  has  the  fora  _u  «  u(x.v)e|,.  The  flow  is  forced  by  a  given  pressure 
gradient  grad  p  -  -Go  .  for  this  situation,  ths  aquation  of  equilibrium  reads 

(2.1)  div(iiVu)  -  -G 
with  the  boundary  condition 

(2.2)  u  ■  0  on  20  . 

Bare  u  is  a  step  function  assuming  the  values  u1  and  u j.  On  the  interface  between  the 
two  fluids  the  velocity  and  shear  stress  have  to  be  continuous.  These  two  conditions  are 
automatically  satisfied,  if  (2.1)  is  taken  in  the  distributional  sense.  These  continuity 
conditions  can  also  be  thought  to  arise  as  natural  boundary  conditions.  The  equations 
(2.1)  and  (2.2)  are  the  Euler  aquations  for  the  variational  problem 

(2.3)  min  /fl  j  «(Vu)2  -  Gtt  . 

uEHq(Q) 

Let  us  denote  the  integral  expression  in  (2.3)  by  r  (u).  Since  this  Is  a  strictly  convex 
functional  of  u,  the  following  lemma  is  immediate. 

Lemma  2. It  For  any  w  CL  (Q)  such  that  |i  >  e  >  0,  there  exist  one  and  only  one 
u  6  Bg(Q),  which  solves  problem  (2.3),  and  hence  problem  (2.1),  (2.2). 

This  lemma  shows  that  the  flow  of  two  fluids  has  a  high  degree  of  nonuniqueness.  It 
says  that  for  any  given  arrangement  of  the  fluids  there  Is  a  corresponding  flow  field. 
Experimentally,  however,  certain  Interface  positions  are  observed,  others  are  not.  The 
question  thus  arises,  which  configuration  of  the  Interface  is  stable. 

McLean  [6]  and  Bverage  [2]  have  suggested  that  the  stable  solution  nay  be  the  one 
which  minimises  (2.3)  also  with  respect  to  |i.  If,  say,  the  volusw  occupied  by  each  fluid 
is  given.  That  Is,  u  is  permitted  to  vary  in  the  class  of  step  functions 
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*  -  {»2  ♦  ,|l1"*l2>xn  *Qic  0  ■M,ur*bl**  lat I  “  ■> 

where  a  la  a  niabir  which  specifies  the  volume  ration  and  xn  denotes  the  character ia- 

a1 

tic  function  of  flji  l.e.,  Xjj  (*>  ■  1  whan  ««lt(  and  xn  "  0  otherwise.  If  we 

denote  by  D  the  rate  of  viscous  dissipation,  and  by  tr  the  Mechanical  work  done  by  the 
fluid,  then 

y«>  ■  i  »-*  • 

tor  any  steady  flow,  we  have  D  -  W,  and  hence  the  flow  that  Mininizes  P  is  the  one 

u 

that  maximizes  D  or  W.  Also  since  W  -  G  u.  It  aaxialsea  /  u,  the  volimm  flux  of 
the  fluid. 

In  general,  we  have  no  proof  that  the  problen 

ain  r  (u) 

1  •* 

V  6* 

always  has  a  solution.  A  nlnlmlzer  does  exist,  however,  if  we  impose  an  a  priori  boundary 

on  the  length  of  the  fluid  Interface.  For  this  purpose,  let  us  define 

*L  -  (u  e  *1  For  any  e  >  0  there  is  a  set  of  at  most  [L/e]  ♦  1 

closed  disks  of  radius  e  which  completely  covers  9(1^}  . 

this  condition  Is  a  mathematical  rlgorlsatlon  of  the  statawent  that  the  length  of  9Q 
la  less  than  or  equal  to  2l. 

Theory  2.2.  (due  to  M.  Rsnardy) 

•  1 

there  exists  u  8  ♦  and  u  6  H.(tl)  such  that 

Xi  0 

* 

F  #(u  )  -  Min  r  (u)  . 

M  U«L  " 

u8h’(0) 

Frooft  bet  tin,  u  be  a  mlnlalslng  sequence.  Then  clearly  {a*}  Is  uniforaly  bounded  in 
and  hence  has  a  weakly  convergent  subsequence.  We  shall  show  next  that  any  sequence 
{|in}  In  *  has  an  almost  everywhere  convergent  subsequence  whose  llMlt  is  again  in  •  • 


To  see  this,  lot  e  >  0  ha  art  van.  For  each  4  ,  va  have  at  Most  lv/c)  *  1  e-disks 


He  can  now 


covering  30*.  Denote  the  centers  of  these  c -disks  by  Hi*  Wo  can  now 

\  \  \  "v 

extract  a  subsequence  H  such  that  ,  Hj  * '  *  * 'N[x./e]^1  convac9*'  ••y«  to 

Hj , . .  •  ,M  jL_/e j •  Then,  for  any  fixed  4  and  k  large  enough,  the  boundary  of  0^  Is 

covered  by  (e+4  ) -disks  centered  at  , . . .  «Mfjye]  +  t  and  thus  by  a  set  of  Measure 
2 

*(c+4)  ((L/C]  ♦  t).  The  complement  of  this  sat  consist  of  finitely  Many  components  on  each 
nk 

of  which  the  u  are  constant  for  k  large  enough.  Hence  we  can  extract  a  convergent 

subsequence  there.  Repeating  this  argument  for  a  sequence  e  such  that  c  ♦  0  and 

v  in 

using  the  standard  diagonal  arginent,  we  find  an  almost  everywhere  convergent  subsequence 

of  {u"}.  Now  assume  u"  e  4,  u"  ♦  u  almost  everywhere  and  u  i  • .  •  Then  there  is  some 

L  '  L 

e  >  0  such  that  the  boundary  of  (u )  *  (<x,y)  e  Q,  u(x,y)  “  u  ^ }  cannot  be  covered  by 

(L/eJ  +  t  e-disks,  if  that  is  so,  then  there  is  a  finite  number  P, ,P2, . . . ,PN  of  points 

in  30^(11)  which  cannot  he  covered  by  tL/e]  ♦  1  e-disks.  However,  for  any  4  >  0,  and 

any  k  “  t,...,N,  there  is  a  point  p*  e  3a"  in  the  4-neighborhood  of  Pk,  if  n  is 

chosen  large  enouqh.  Since  u"  e  ♦  ,  it  follows  tKat  the  p"  can  be  covered  by 

L  K 

ll/t]  ♦  1  e-disks,  and  hence  that  Py  are  covered  by  [L/e]  +  1  e+6-disks.  Letting 
4*0,  this  yields  a  contradiction,  thus  by  extracting  convergence  subsequences,  »e  find 
U*  ♦  |i  a.e.  and  u"  ♦  u  weakly  in  H^(tl ) .  It  is  an  easy  consequence  that 
nn  »  (un)  1  *  n  *  |i  1  .  Moreover  there  is  a  sequence  vn  of  convex  comblnatins  of  the 
un  such  that  v°  *  u  strongly  in  H^IQ)  and  Tv"  ♦  Vu  almost  everywhere.  Since  the 


functional 


is  convex,  we  have 


K(n,u)  -  /g  n  ^(Vu)2  -  Gu 


lie  K(nn,vn)  <  lim  K(nn,un)  , 


where  n"  is  an  appropriate  convex  combination  of  the  n”.  According  to  Patou's  lemma,  we 


F  (u)  «  K(n,u)  <  11m  k(on,vn)  <  lim  K(nn,un)  ”  lim  F  (un) 
‘  n+-  n*»  u" 


This  proves  that  the  pair  (ji,u)  is  a  minlxiiear. 
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For  the  apecial  com  of  a  circular  pipe,  an  eleuentary  argument  la  poaalbla.  The 

following  nodi  flea  an  idea  of  Kverege  (2].  Lat  Q  be  a  circular  disk,  and  lot  Ug  bo  tho 

**0  — 

solution  of  AUg  “  -G,  Ug|Ja  *  0.  Ha  put  u  *  ~  ♦  u.  Than  tho  integral  in  (2.3)  la  equal 
to 

fa  Vu0Vu  -  GU 

♦  /Qj  5  ^(VU)2  -  ^  (7u0)2 

♦/QiiM,<Vu)2-^.(7u0)2 

wharo  flj  -  Q\R1  la  the  region  occupied  by  high  viscosity  fluid,  the  first  of  the 
expressions  la  aero.  The  ten 

/«2  2i72  <Vuo)2  +  4,  ^  (V»o>2  *• 

2 

aaxleal  if  Q^,  the  region  of  higher  viscosity,  is  where  (VuQ )  takas  lta  sealleat 

values >  i.e.  If  Qj  la  a  disk  in  the  center  of  the  pipe.  In  this  case  301  la  a  line  on 

which  Ug  la  constant.  If  we  then  choose  u  «  0  in  the  outer  region,  and  u  »  const,  in 

the  inner  region  such  that  u  is  continuous,  then  the  continuity  of  velocity  and  shear 

stress  across  the  Interface  is  satisfied.  The  expression 

\  J«2,V»)2+/ni?'V’»>2 

becomes  sero,  which  is  clearly  its  einlnal  value.  He  have  therefore  proved 

Theorem  2.3. 

If  8  is  a  circular  disk,  then 

■in  /  j  u(?u)2  -  Gu 

umj(a> 

we* 

is  attained  for  a  concentric  configuration  with  the  higher  vlacoaity  in  the  inner  region. 
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3.  Lin— r  Stability 

W»  consider  •  circular  plpa  of  radius  Rj  and  a  baalc  flow  given  by  an  axial  velocity 


<«-m>R2  -  r2]  0  <  r  <  »1 

*2  “  4^  <R2  _  t2)  R1  <  r  <  *2 


Hara  u 1 »  u2  ara  the  vlacoaitiea  of  the  fluid*  in  the  inner  and  outer  region,  and  G  ia 
the  applied  pressure  gradient.  He  auperlapose  an  inf initeslmal  diaturbance  which  la 
periodic  In  the  axial  direction  and  uae  a  Fourier  expanaion  In  the  azimuthal  direction)  we 
therefore  represent  the  diaturbance  by 

(u,v,wjp)exp(-lact  ♦  iaz  +  in  8) 

where  u,  v,  w  denote  the  radial,  azimuthal  and  axial  velocities  and  p  is  the 
pressure,  the  stability  equations  are  derived  from  the  Navier-Stokes  equations  plus  the 

✓ 

conditions  that  velocities  and  tractions  are  continuous  at  the  interface.  From  these 
conditions  the  following  eigenvalue  problem  for  c  is  obtained  [4J  (for  simplicity,  we  put 
the  density  of  both  fluids  equal  to  1): 

2 

ia<lf-c)u  -  -p'  *  (ruT  -  ~~  u  -  o2u  - 

r  r 


(3.2) 


lo(W-c)v  -  -  ~E  u  (i  (rv')*  -  S-j!  V  -  02v  +  ^iSS] 

r  r 

1  2 

io(W-c)w  ♦  uW’  ■  -lop  +  g [—  (rw1 )'  -  ~  w  -  o  w] 

r 


(ru)'  ♦  inv  +  riaw  «  0 


The  boundary  and  interface  conditions  are 


-8- 


r 


u  -  v  “  w  -  0  at  r"*2 
tvl  ■  lul  -  0  'l 


(3.3) 


[p]  -  2biu‘] 

CM v* J  ♦  (lnu-v)lii)  •  0 
*1 

Iu«’]  ♦  lau(w)  “  0 

ulw* )  +  (M-c)la {«]  -  0  . 


at  r  -  R. 


Hara  a  prlae  danotaa  jj,  1  ]  danotas  tha  dlffaranca  between  tha  values  of  a  quantity  on 


both  aidaa  of  tha  interface,  and  M  la  given  by  (3.1). 


4.  Numerical  calculations 
4.1.  Method  of  calculation 

Me  ellainata  w  and  p  froa  tha  third  and  fourth  aquation  in  (3.2),  so  that  tha 

problaa  reduces  to  finding  u  and  v  froa  tha  rasuilning  aquations  and  tha  boundary 

conditions  (3.3).  Mhen  eliminating  w  froa  tha  fourth  aquation,  we  divide  by  r,  thus, 
even  if  we  Inpose  the  condition  that  u  and  v  are  saooth  at  r  »  0,  we  allow  w  to 

have  an  (unphyalcal)  ^  -  singularity.  Hie  third  aquation  in  (3.2)  has  a  regular  singular 

point  at  0,  and  -1  is  a  root  of  tha  lndiclal  aquation  precisely  for  n  “  1.  Since  the 
second  root  la  >1,  tha  unbounded  solution  generally  behaves  like  —  ♦  ar  In  r  near 
r  ■  0.  it  can  be  shown  that  tha  coefficient  a  of  tha  logarithmic  tern  vanishes  iff 
c  "  Wj(0)  -  3n jdl .  This  appears  in  our  calculations  as  an  extra  eigenvalue  for  n  «  1, 
which  aust  be  dlsaissad  aa  unphyalcal. 

For  siaplicity,  we  take  the  inner  radius  R^  equal  to  1.  Me  expreas  the  radial 
dependence  of  the  disturbance  as  a  linear  combination  of  Chebyshev  polynomials  TB(r)  - 
cos(a  arccos  r)  aa  suggested  by  Orssag  (8)  and  truncate.  Hence  the  radial  velocity  is 

H 

u  -  l  u^T^tr)  0  <  r  t  1 

B“0 

N 

U  •  l  U  T  (r)  1  4  r  4  s 

- a  «I  Bl  « 


9- 


where  r  -  -1  ♦  rufM  from  -1  to  1.  Bio  torso  of  highest  differential  order 

*2  1 

ere  r4u*  •••  and  r3v*  •  *  In  the  flret  equation  of  (3.2)  end  r3u* "  and  rSv**  In  the 
eecond  equation.  We  would  like  the  terse  end  to  be  approx lse ted  to  the  ease 

degree.  Since  we  choae  a  to  be  an  Nth  degree  polynomial,  both  terse  are  of  degree 
N  -  4  if  v  la  a  polynomial  of  degree  N  -  1.  Ne  thus  pot 

N-1 

v  -  I  v  T  (r)  0  <  r  <  1 

~o  ■■ 

N-1  _ 

v  -  l  v  T  (r)  1  «  r  «  R„  . 

J;0  ■■  2 

Zf  the  asisothal  number  n  la  even,  then  u  and  v  moat  be  odd  functions  of  r 
(coneoqoently  p  and  w  are  even),  if  n  la  odd,  then  u  and  »  are  even.  Coneletent 
with  thle,  the  eons  In  the  inner  region  are  restricted  to  odd  (or,  rasp.,  even)  Chebyehev 
polynomials,  the  expressions  for  a  end  v  are  inserted  into  the  equations,  which  are 
then  truncated  at  orders  N  -  4,  N  -  3  rasp.  For  convenience,  we  take  N  to  be  odd.  In 
the  outer  region,  we  then  hews  2N  ♦  1  unknowns  and  2N  -  5  equations.  In  the  Inner 
region,  we  get  N  ♦  1  unknowns  and  N  -  2  equations  If  n  Is  odd,  while  we  get  N 
unknowns  and  N  -  3  equations  If  n  is  even,  together  with  the  9  boundary  and 
Interface  conditions,  this  yields  3N  ♦  2  (for  n  odd)  or,  Tf*  ♦  1  (for  n  even) 
equations  and  unknowns. 

The  eigenvalues  c  of  the  resulting  matrix  equation  were  computed  in  complex  double 
precision  on  a  VAX/VN8  V02  system  using  the  INK.  routine  BIQXC.  the  flow  Is  unstable  if 
the  Imaginary  part  of  c  is  positive. 


4.2.  accuracy  and  convergence 

Ne  compared  the  eigenvalues  for  the  one-fluid  case  -  u2  with  the  results  of 
Salwen  and  Groach  [9]  and  Salwen,  Cotton  and  Groech  [10].  For  n-1,  they  list  four 


eigenvalues  with  the  smallest  Imaginary  parts  at  Reynolds  number 


-  100 


and 
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higher,  to  5  deelaal  place#.  the  velocity  ecale  -j~  la  fixed  at  1,  aad  OH^  •  1.  We 
«•  Mainly  lntaraatad  la  low  Haynolda  ntafeera,  ao  wo  ehaokad  our  oaloulatlona  at  Mo  -  100 
trunoated  at  M  -  1»  for  varloua  value a  of  1]  to  bo  good  la  blooaponant  flown,  A 
conparleoa  la  abowa  la  figure  1.  for  other  aalaatbal  aodaa  aad  lower  valuaa  of  Me,  our 
eigenvalue a  ware  ohaokad  agalnat  graph a  publlahad  by  falwan  at  al.  At  beat,  we  oaa  read 
their  graphe  to  2  digit#,  aad  our  raaulta  for  W  ■  it  agree  within  that  accuracy.  Our 
valuaa  for  He  «  10  aad  100  are  ahowa  la  figure#  2  aad  2. 


figure  1 

He  -  QH^/4ua  »  100  , 

Velocity  aealo  -  3R^/4y  -  i 
node  1 


*2 

D 

s 

.04 

elgenvaluea 

Publlahad  by  ^ 

Halwen,  Cotton  a  Oroaohj 

1 

.01 

.57256, -.14714 

.55198, -.37446 

.78735, -.47946 

.66247, -.74907 

Our  oonputatlon 

El 

EM 

g 

a 

m 

■ 

.66248, -.74908 

with  H  -  19,  -  Uj 

I"-. 

.02 

m 

.55 199,-. 37446 

■ 

.66248, -.74907 

il  .28 

1 _ 

B 

.32 

.57254,-. 14716 

.55214, -.37448 

.78736, -.47950 

.66220, -.74907 

Coaiparleon  of  our  one-fluid  tu1  -  w2>  elgenvaluea  with  thoee  publlahad 
by  Halwen  at  al  for  node  1  and  Ha  -  100. 


rlgura  2 


*2  -  1.0 

»1  ■  U2  ■  0.16 

Rj  -  1.6 

G  -  0.24 
N  -  19 


Mode 

algaaraluoa 

0 

.66676,-1.5776 

.81338,-  2.6913 

.66657,-  5.0195 

.72389,-  7.1393 

.66666,-10.448 

1 

.49106,-1.3935 

.76202,-  2.8073 

.62772,-  4.8498 

.70556,-  7.2282 

2 

.40791,-2.6274 

.76307,-  4.2327 

.57022,-  7.0840 

.71238,-  9.6639 

3 

.36440,-4.0905 

.73178,-  5.9084 

.53554,-  9.5574 

.70101,-12.372 

4 

.33513,-5.7962 

.69400,-  7.8343 

.50882,-12.289 

.68325,-15.348 

5 

.31324,-7.7438 

.65689,-10.006 

.48652,-15.281 

.66345,-18.587 

Oar  on«-(liid  aiganvaluai  for  aodaa  0  to  5  at  Ra  ■  10  which  can  ba 
coaparad  with  publlahad  grapha  of  Salwon  at  al 
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rigor*  3 


Da  -  100  a*j  -  1.0  W19 


1.33 

• 

2.0 

3.0 


.0123 

.02 

.05 


1.251.0123 

•  i 

2.0  .02 
5.0  j .05 
- 1 - 


1.25  .0123 


2.0  .02 


5.0  !.0S 


1.25  .0125 

•  i 

2.0  |  .02 

5.0  !.05 


1.25  1.0125 


2.0 

5.0 


.02 

.05 


.032 

.02 

.000 


.032 

.02 

.000 


.032 

.02 

.008 


algaawaloas 


.68626, -.27345 
.88627,-. 27345 


.49400, -.25600 

.49400, -.25600 


.032 

.02 

.008 


.032 

.02 

.008 


.44087, -.39187 
.44088, -.39189 


.40024,-. 55701 

.40024,-. 55701 
.40023, -.55702 


.36772, -.74937 

,36772, -.74934 


81160,-  .28413 

,81160,-  .28413 


,70711,-  .62217 


.70710,- 

.70710,- 


.62217 

.62218 


.64116,-  .71996 


.64097,. 

.64096,- 


.71991 

.71990 


.62004,-  .88366 


.62006,- 

.62007,- 


.88365 

.88361 


.59648,-1.0844 


.59647,- 

.59644,- 


1.0845 

1.0845 


.64648, 

-  .45939 

.64647, 

-  .45937 


.56483, 

-  .58500 

.56482, 

-  .58498 
.56482, 

-  .58497 


.56903, 

-  .89443 

.56913, 

-  .89404 
.56914, 

-  .89402 


.53298, 

-1.1892 

.53294, 

-1.1892 

.53301, 

-1.1892 


.72999, 

•  .73274 

.72993, 

-  .73276 


.66679, 

-1.0274 

.66663, 

-1.0277 


.68056, 

-1.0059 

.68051, 

-1.0061 

.68052, 

-1.0062 


.67818, 

-1.2737 

.67734, 

-1.2767 

.67722, 

-1.2769 


.50485, 

-1.4990 

.50473, 

-1.4986 

.50458, 

-1.4987 


.66487, 

-1.5714 

.66482, 

-1.5714 

.66372, 

-1.5719 


.64873, 

-1.8911 

.64868, 

-1.8919 

.65052, 

-1.8920 


.69394, 

-1.3568 

.69444, 

-1.3579 


Oar  on*-fluld  alganvaloaa  for  aodaa  0  and  2  to  5  at  Ra  ■  100  which 
can  b*  cowparad  with  publlahad  graph*  of  Sal wan  at  al. 


-13- 


Dm  Instability  dan  to  a  junp  In  viscosity  la  governed  by  an  lnterfaeial  node,  which 


at  »1  -  tij»  la  c  •  WOt,),  aa  can  ba  aaan  fron  tha  last  lntarfaoa  condition  In  (3.3). 
tha  value  of  e,  and  In  particular  Its  Imaginary  part,  la  typically  aaallar  by  ordara  of 
nagnitude  than  tha  resMining  eigenvalues.  Bence  e  will  generally  ba  approx iaa tad  laas 
wall  than  tha  othar  eigenvalues.  Since  wo  ara  aainly  Intaraatad  in  tha  alga  of  In  o,  wa 
ara  aatlaflad  with  2-digit- accuracy. 

Appropriate  Noynolda  nuabara  for  tha  innar  and  outer  fluids  ara  b  t  • 
and  Rs  2  -  *2 < R1 > R2^ 2  raapaetlwaly.  In  oar  figures,  ha  will  danota  a  rafaranoa 
Reynolds  nunber  wj(°>*2^i2  in  with  tha  ona-flald  flow,  and  'velocity  scale"  will 

rafar  to  W3(0).  tha  into)  la  proportional  to  tha  valoelty  acala  so  wa  have  ehoaan 

Wj(0)  -  1. 

Figure  4a  to  4a  display  saaplas  of  convarganca  testa,  thara  ara  three  aaln 
features.  First,  since  wa  used  an  equal  n unbar  of  nodes  in  tha  Inner  and  outer  regions,  we 
expect  convergence  rates  to  worsen  as  1,/lj  novas  away  fron  O.S.  Figure  4a  oenparas 
*l/«2  “  0.5  with  0.2.  At  R1/R2  “  O.S,  all  awdas  have  converged  to  at  least  2  digits 
by  K  -  19.  At  R1/R2  “  0.2,  *  -  39  la  required  for  slnllar  accuracy  for  nodes  t  to 
4.  Also,  as  l^/Fj  *  0,  in(c)  ♦  0  so  round-off  errors  enter.  Our  nethod  of  expansion  has 
an  accuracy  of  infinite  order  so  that  once  lu(c)  has  converged  to  2  digits,  convergence 

to  aore  digits  will  occur  for  values  of  H  which  ara  not  such  larger. 

Secondly,  at  low  Re,  convergence  is  worse  for  the  higher  Bodes  because  la(c) 

decreases  fast  In  nagnitude  as  node  n unbar  increases.  Refer  to  figure  4b  for  Rj/Rj  “  0.6, 

oRj  >0.1,  Re  >  0.1.  However,  computations  need  not  be  done  for  very  snail  Re  because 
ln(c)  is  proportional  to  Rs  In  the  Unit  as  Re  ♦  0.  Figure  5  shows  that  this  low  Rs 
ragine  extends  up  to  about  Re  1,  Re  2  -  10,  or  Re  -  1  for  the  values  of  we 

consider.  For  higher  Re  (figure  4c)  the  nagnitude  of  in(c)  is  similar  for  all  nodes. 
Rota  that  Rs  1  Increases  as  li  ^/U j  decreases  and  this  accounts  for  tha  slower 
convarganca  at  snaller  U ,/Mj/  (see  figure  4d).  At  Re  ■  1000,  the  nagnitude  of  tha 
lntarfaelal  node  la  slnllar  to  that  of  the  stable  eigenvalues  associated  with  one-fluid 
flow.  Convergence  lnproves  for  given  a  and  R,/R2  ••  Increases. 
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Cmwrvmo*  Uiti 
Velocity  Seal*  ■  1 

Flgora  4a 

Ha  -  100  ,  «Rj  “  1 


.1  |3.3199,  -.809088*1 
1.0454,  -.994401-1 
.75041, -.283951-3 
.74972,  .274421-3 
1.1042  ,-.321S0C-2 
.94904, -.245591-3 
.94209,-. 14S941-3 
.95824,  .143438-3 


2.3200,  -.809408-1 
1.0494,  -.994598-1 


.94234,  .4459888-9 
.95804,  .441098-4 


1.0021  ,-.931418-3  (8-25)  .94519,  .199518-3 
1.0432,  .801448-2 


.75041, -.2838 58-3 
.74972,  .276748-3 
1.1062  ,-.321658-2 
.94904, -.245418-3 
.96236,  .161578*4 
.95805,  .3487748-4 
.95995,  .112748-3 


;  i 

! 

I 

■  i 

Flgura  4b 

Ha  -  0.1  ,  a*j  ■  0.1 

R,/!^  aoda  u  1/U2 

o 

|  1 

8-19 

8-23 

0.8  0 


.93533, -.21227  1.0734,  -.22125 

.47541, -.944398-1  .47743, -.953788-1 

.44044, -.282488-1  .44029, -.281758-1 


.46588,  .610558-1 
.44054,  .392698-1 
.42812,  .207658-1 


.46955,  .514548-1 
.43580,  .365658-1 
.41937,  .263938-1 


8,/Bj  Loda 


Flyura  4d 
100  a8-  »  1 


«<**■,  *^T— <tr-' 


rigor*  4* 

Ra  -  100  aRj  -  10 


R,/Rj 

aod* 

v1A>2 

e 

W19 

8-23 

0.6 

0 

.2 

.74558, .195351-1 

.74624,. 192 168-1 

.4 

.67284, .866058-2 

.67283,. 868108-2 

1.2 

.63835, .373678-3 

.63835, .373678-3 

4 

.2 

.71083,. 280608-1 

.71 152,. 280258-1 

.4 

.65854, .946028-2 

.65853, .946718-1 

1.2 

.63931,. 403208-3 

.63931, .403208-3 

0.9 

0 

.2 

.30635, .232488-1 

.30454,. 2095 78-1 

.23388,. 109088-1 

.23385, .107588-1 

1 

.30903, .232558-1 

.30379,. 214458-1 

n*  «  1 


R1/R2 

aod* 

“1*2 

c 

H-19 

8-23 

0.8 

0 

.2 

.38505, .103638-1 

aaa*  to  all  digit* 

Conrargance  t**t*  for  velocity  seal*  1,  R,  -  1,  c  1*  th* 
intarfacial  aiganvaloa,  R  ia  dafinad  in  {4. 
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rlfwt  s 


Tibia  of  o 


Mod* U/Uj  Ra  -  0.1 


oRj  "  1.0 

R  -  1.6 

If  U,/V2  -  1  ,  c  -  .64 

Ra  -  1.0  1  Ra  ■  10  I  Ra  -  100 


101101-2 

4S676I-4 

263741-5 

544301-6 

785881-6 

140191-5 

184281-5 

215181-5 

236701-5 


,  .173131-2 
,  .303671-3 
.  .945151-4 
.  .291091-4 
,-.156651-4 
,-.250331-4 
,-.309721-4 
,-.348591-4 
,-.374331-4 


.97610, 

.80554, 

.72215, 

.67266, 

.61692, 

.59981, 

.58667, 

.57630, 

.56794, 


.989761-2 

.455701-3 

.263301-4 

-.544651-5 

.785811-5 

.140191-4 

.184281-4 

.215181-4 

.236701-4 


.172721-1 
.303621-2 
.945121-3 
.291121-3 
-. 156651-3 
-.250311-3 
-.309711-3 
-.348581-3 
-.374331-3 


1.0332, 

.80832, 

.72243, 

.67270,* 

.61691, 

.59980, 

.58666, 

.57628, 

.56793, 


.275461-1 

.362051-2 

.220821-3 

.579551-4 

.789101-4 

.140471-3 

.184481-3 

.215311-3 

.236781-3 


297831-1 

939091-2 

289551-2 

156021-2 

249491-2 

308911-2 

347881-2 

373741-2 


.98511,  . 
.74489,  . 
.67945,  . 
.65139,  . 
.63404,-. 
.62850,-. 
.62287,-. 
.61731,-. 
.61193,-. 


21492 

442671-1 

944791-2 

227101-2 

103151-2 

162381-2 

200421-2 

226021-2 

243531-2 

23978 

11063 

555211-1 

216041-1 

115361-1 

180111-1 

223171-1 

255281-1 

280671-1 


2  .2  .70888,  .606571-3  .70902,  .606251-2  i  .72345,  .557101-1  ; .84503,  .13246 

.4  .68357,  .101091-3  .68358,  .101081-2  .68385,  .100501-1  1.70348,  .568981-1 

.6  .66499,  .299531-4  .66499,  .299521-3  .66493,  .298441-2  ! .66489,  .230131-1 

.8  .65092,  .876281-5  .65092,  .876141-4  .65089,  .873171-3  ' .64945,  .704491-2 

1.2  .63138, -.423081-5  .63138, -.423481-4  .63139, -.422221-3  .63257, -.3380 71-2 

1.4  .62447, -.640901-5  .62447,-. 641191-4  .62449, -.639401-3  j .62638, -.513371-2 

1.6  .61886, -.7521 11-5  .61886, -.751911-4  .6 1890, -.749931-3  ] .621 17, -.6051 51-2 

1.8  .61428, -.802211-5  .61428, -.802201-4  .61432, -.800181-3  .61677,-. 649201-2 

___  2.  .61051, -.816591-5  .61051, -.816511-4  .61054,-. 814581-3  .61305,-. 664101-2 

3  .2  .68296,  .271771-3  .68300,  .271691-2  .68732,  .262461-1  F. 76245,  .897821-1 

.4  .66619,  .399521-4  .66620,  .399531-3  .66629,  .398361-2  .67516,  .294101-1 

.6  .65453,  .993841-5  .65453,  .993801-4  .65452,  .992121-3  .65474,  .866951-2 

.8  .64615,  .238721-5  .64615,  .238721-4  .64615,  .238301-3  .64595,  .209131-2 

1.2  .63540, -.731841-6  .63540, -.731791-5  .63S40,-. 730161-4  .63551, -.6 12391-3 

1.4  . 63191, -.831711-6  .63191, -.831021-5  ,  .62925, -.665521-4  .63205,-. 666361-3 

1.6  .62924,-. 668151-6  .62924, -.667721-5  .62925, -.665521-4  .62935, -.492641-3 

1.8  .62719, -.397631-6  .62719, -.397721-5  .62719, -.395211-4  .62726, -.227351-3 

2.0  .62S6t, -.897121-7  .62561, -.898921-6  .62 561, -.877851-5  .62562,  .681791-4 
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.2 

— 

.66341, 

.136838-3 

.66343, 

.136818-2 

.66490, 

.134448-1 

.4 

.65379, 

.173348-4 

.65379, 

.173338-3 

.65303, 

.173078-2 

.6 

.64742, 

.326188-5 

.64742, 

.324588-4 

.64742, 

.324308-3 

.8 

.64306, 

.459848-6 

.64306, 

.455708-5 

.64306, 

.455308-4 

1.2 

.63781, 

.152728-6 

.63781, 

.152308-5 

.63781, 

.152108-4 

1.4 

.63623, 

.470258-6 

.63623, 

.458168-5 

.63623, 

.457958-4 

MEM 

.63907, 

.782298-6 

.63507, 

.779238-5 

.63507, 

.778968-4 

■  I 

.63422, 

.106238-5 

.63422, 

.107198-4 

.63422, 

.107178-3 

I  1 

.63360, 

.132548-5 

.63360, 

.132508-4 

.63359, 

.132468-3 

.2 

.65156, 

.740448-4 

.65157, 

.740378-3 

.65219, 

.733388-2 

.4 

.64662, 

.812208-5 

.64662, 

.812218-4 

.64664, 

.811558-3 

.6 

.64348, 

.996728-6 

.64348, 

.997588-5 

.64349, 

.997528-4 

.8 

.64141, -.615008-7 

.64141, -.622228-6 

.64141, -.621508-5 

1.2 

.63902, 

.283918-6 

.63902, 

.282948-5 

.63902, 

.282848-4 

1.4 

.63833, 

.573938-6 

.63833, 

.575428-5 

.63833, 

.575208-4 

1.6 

.63705, 

.827568-6 

.63785, 

.827098-5 

.63784, 

.826688-4 

1.8 

.63750, 

.103288-5 

.63750, 

.103058-4 

.63750, 

.103028-3 

2.0 

.63725, 

.118778-5 

.63725, 

.119078-4 

.63725, 

.119048-3 

.71072, 

.65797, 

.64762, 

.64314, 

.63774, 

.63609, 

.63487, 

.63398, 

.63332, 


.67915, 

.64867, 

.64381, 

.64151, 

.63895, 

.63821, 

.63769, 

.63731, 

.63705, 


.625088-1 

.148038-1 

.305188-2 

.430598-3 

.150038-3 

.447468-3 

.759688-3 

.104508-2 

.129238-2 


.434358-1 

.747958-2 

.991028-3 

.505758-4 

.270928-3 

.554268-3 

.799578-3 

.999138-3 

.115708-2 


Tab  la  of  intarfacial  aiganaaloaa  e  at  aodirat*  oRJf  allowing  tha  llnaar 
bahaviour  of  a  for  aaall  8a. 


rigor*  6 


aR. 


0.1 


VaX.  seal*  -  1 

R  -  1.6 

If  "  1 


c  -  .64 


u1/U2 

Ra  -  0.1 

»•  -  1.0 

Ra  -  100 

M3 

N-21 

9"  19 

.2 

1.0278,  .768321-4 

1.0279,  .768138-3 

1.0789,  .115378-1 

.4 

.82497,  . 1 5796R-5 

.82497,  .157948-4 

.82722,  .833758-3 

.73040, -.343 198-6 

.73040, -.343 148-5 

.73063, -.377598-3 

.67567,-. 182628-6 

.67571, -. 1856oC-3 

.61490,  .116028-6 

.61490,  .115998-5 

.61489,  .116308-3 

1.4 

.59629,  .186668-6 

.59629,  .186768-5 

.59628,  .187038-3 

MKm 

.58193,  .230248-6 

.58193,  .230568-5 

.58192,  .230768-3 

mm 

.57052,  .258578-6 

.57052,  .258198-5 

.57051,  .258338-3 

mm 

.56123,  .275788-6 

.56123,  .275848-5 

.56122,  .275948-3 

.2 

.73041,  .138808-3 

.73041,  .138808-2 

.76520,  .12215 

.4 

.69982,  .283S5B-4 

.69982,  .283528-3 

.69829,  .281838-1 

.6 

.67567,  .996968-5 

.67567,  .996668-4 

.67466,  .981138-2 

.8 

.65614,  .332408-5 

.65614,  .332138-4 

.65578,  .327208-2 

1.2 

.62645,-. 193528-5 

.62645,-. 193518-4 

.62664,-. 191648-2 

1.4 

.61491, -.315598-5 

.61491, -.315478-4 

.61519, -.313108-2 

1.6 

.60496, -.395S4S-5 

.60496, -.395678-4 

.60529, -.393378-2 

1.8 

.59629, -.449658-5 

.59629, -.449628-4 

.59664, -.447598-2 

uL _ 

.58868, -.485948-5 

.58868, -.486248-4 

.58903, -.484538-2 

.2 

.71491,  .596748-4 

.71491,  .596778-3 

.72743,  .561258-1 

.4 

.68751,  .105528-4 

.68751,  .105518-3 

.68756,  .104948-1 

.6 

.66732,  .331228-5 

.66732,  .331198-4 

.66720,  .329558-2 

.8 

.65196,  .101258-5 

.65196,  .101288-4 

.65192,  .100798-2 

1.2 

.63050, -.518178-6 

.63050, -.519528-5 

.63053, -.517568-3 

1.4 

.62285, -.803428-6 

.62285,-. 804418-5 

.62289, -.801548-3 

1.6 

.61661, -.962448-6 

.61661, -.961878-5 

.61666, -.9590 38-3 

1.8 

.61 148,-. 104558-5 

.61148,-. 104598-4 

.61153,-. 104268-2 

2. 

.60722,-. 108328-5 

.60722,-. 108348-4 

.60727, -.1080 58-2 

.2 

.68607,  .273408-4 

.68607,  .273468-3 

.69023,  .265128-1 

.4 

.66817,  .416538-5 

.66817,  .413748-4 

.66825,  .412658-2 

.6 

.65566,  .899128-6 

.65566,  .107668-4 

.65565,  .107028-2 

.8 

.64665,  .328508-6 

.64665,  .269928-5 

.64664,  .269458-3 

1.2 

.63501, -.989268-7 

.63501, -.947258-6 

.63501, -.933318-4 

1.4 

.63121, -.103268-6 

.63121,-. 115668-5 

.63121,-. 116308-3 

1.6 

.62829, -.170038-6 

.62829,-. 110258-5 

.62829,-. 108358-3 

1.8 

.62603, -.279778-6 

.62603, -.840668-6 

.62604, -.864948-4 

2. 

.62428, -.490648-7 

.62428, -.518888-6 

.62429, -.587438-4 

-20- 
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.2  .66492, 

.4  .65470, 

.6  .64792, 

.8  .64327, 

1.2  .63765, 

1.4  .63595, 

i  1.6  |  .63470, 

I  1.8  !  .63379, 

j  2.  j  .63311, 


.138241-4  .66492, 
.177628-5  .65470, 
.347298-6  .64792, 
.520898-7  .64327, 
.109658-7  .63766, 
.437858-7  .63595, 
.684508-7  .63470, 
.988728-7  .63379, 
.133598-6  .63311, 


.138228-3  .66646, 
.178058-4  .65474, 
.344758-5  .64792, 
.527008-6  .64327, 
.112578-6  .63765, 
.398038-6  .63595, 
.718588-6  .63470, 
.974518-6  .63378, 
.125668-5  .63311, 


.135968-1 

.177598-2 

.344138-3 

.526218-4 

.112168-4 

.396728-4 

.706988-4 

.993738-4 

.124548-3 


table  of  lntarfaclal  eigenvalues  at  aoderate  Re,  showing  the  linear 
behaviour  of  c  for  ssutll  or,  and  a  Re. 


Bioko* 


thirdly,  wo  boro  chocked  with  Blckox's  figure  2  (1971)  whoro  a  lo  swell . 
o bowed  thet  given  p  ^Aij  end  R^/Rj  and  snail  enough  a,  results  for  any  Me  can  bo 
deduced  f roe  calculations  et  one  value  of  Met  nenely  thet  in(e)  la  proportional  to 
a  s  Reynolds  nunber  and  real  (o)  le  not.  Hie  figure  2  shove  that  ln(c)  versus 
Mj/Mj  <  1  for  selected  values  of  p  <  1  for  nodes  0  and  1.  the  qualitative 
behaviour  In  figure  6  egress  with  Blckox's  graph.  Direct  ocnparlson  can  be  node  at 
p “  0.2.  In  hie  notation,  our  ln(e)  is  atc^/Dw^O).  Be  conputee  (c^/lR)  *  104 
where  R  “  Wj.'OIRj/Vj  so  that  with  our  values  at  Re  ■  0.1,  ”  °*3>  R^/Rg  “  0.6  In 

figure  6,  we  obtain  (c^/iR)  *  104  ■  7.17  for  node  0  and  13  for  sode  1  which  agree 
with  hie  graph,  the  deviation  fra*  linearity  in  Re  at  snail  p^Ai^  end  Me  -  100  le 
due  to  the  fact  that  Rlckox'a  result  la  a  "eeall  a  x  Reynolds  nwdser*  result,  thus  hie 
application  to  oil-water  flows  when  the  Reynolds  rnstber  Is  of  order  10s  would  require 
extraattly  snail  a . 

For  large  aR^,  the  response  le  very  einiler  for  ell  nodes  so  that  convergence  tests 
are  done  for  one  or  two  nodes.  See  figure  4e.  The  ln(c)  is  largest  at  P1/P2  «.2  where 
convergence  is  worst  for  our  range. 

5.  Results 

When  the  Reynolds  nunber  la  sero,  we  can  neglect  the  left  hand  side  of  the  fivst  <hree 
equations  in  (3.2).  Then  c  occurs  only  In  the  Interface  condition  (3.3).  Sy 
substituting  v  -  lv,  w  -  lw,  ws  can  naks  the  equations  (3.2),  (3.3)  real,  and  since  c 
Is  a  sinple  eigenvalue.  It  nust  be  reel.  Thus  I*  c  vanishes  as  Re  ♦  0.  Moreover,  the 
equations  are  Invariant  under  the  change  P  ♦  Xp,  p  ♦  Xp,  c  ♦  Xc,  We  Xw  for  any  positive 
factor  X.  Therefore,  If  the  velocity  scale  and  p  are  changed  by  the  sane  factor  to 
keep  Ms  fixed,  c  Is  proportional  to  the  velocity  scale.  At  low  Reynolds  nisaber, 
therefore,  Ms  c  Is  proportional  to  the  velocity  scale,  while  In  c  is  proportional  to 
the  velocity  scale  tines  the  Reynolds  nwber.  Since  we  have  nornallsed  the  velocity  scale, 
we  have  I*  c  -  Re,  while  Rs  c  tends  to  a  finite  Unit  as  Re  *  0.  This  behaviour  of 
the  Interfacial  node  Is  very  different  fro*  that  of  the  renal ning  eigenvalues,  which  are 


r 


proportional  to  •jjj  aa  Ha  ♦  0.  At  vary  low  Reynolds  nuabera,  thla  laada  to  aavara 
numerical  problaaai  however ,  it  aaaaa  that  wa  can  make  Ha  aaall  enough  to  be  in 
the  "linearised”  rang#  (cf.  chapter  4). 

On  all  grapha/  dark  pointa  rapraaant  confuted  valuaa  and  llnea  which  join  than  are 
interpolants.  Solid  linaa  are  eoaatimes  uaad  to  avoid  confuaion  but  they  are  aleo 
interpolanta.  Humber  a  printed  naxt.  to  the  curve  a  denote  the  asimuthal  mode  number  except 

w, 

on  grapha  7  to  11  where  they  denote  the  ratio  — .  In  all  grapha,  the  velocity  acale  la 

u2 

unity,  on  all  grapha,  except  for  grapha  12  and  25,  the  vertical  axia  meaaurea  the 
imaginary  part  of  the  intarfaclal  eigenvalue  c. 


Moderate  aR^  (~  1),  moderate  Ha  (~  100),  *  [0.2,2] 

The  aign  of  in  (c)  at  given  R^Rj,  doc*  not  *PP**r  to  be  very  aenaitlve  to 

chan  gee  in  aR^  or  Ha.  for  example,  at  Rj  “  1.6,  coapere  Ha  -  1,  a  »  0.06  (graph  1) 
Ha  ■  1,  a  »  0.6  (graph  2),  Ha  -  10,  a  -  0.6  (graph  3)  and  Re  ”  100,  a  -  0.6  (graph  4). 

At  *2  “  1*25»'a  "  0.8,  compare  Ha  ■  1  and  100  (grapha  5,  6). 

The  crucial  paraaetera  determining  the  aign  of  im(c)  are  R1/R2  and  u 1 /u 2 . 

Grapha  7  to  11  dlaplay  the  behaviour  of  ia  <c>  for  a odea  0  to  4  at  Re  ■  100 

a  function  of  R^/Rj,  keeping  —  conatant.  Rote  that  In  many  caaea,  nodea  0  to  4 


aa 


are  not  alaultaneoualy  atabla  or  unatable.  Zm  (c)  decreaaea  to  0  aa  R^Rj  decreaaea 

M1 

from  0.5  to  0  or  increaaaa  from  0.9  to  1.0.  On  graph  7,  the  curve  for  —  “  0.2  ia 

plotted  froa  tracking  the  lnterfadal  eigenvalue  from  R,/Rj  «  0.9.  There  ia  a  degeneracy 

M  g 

for  mode  0  between  R^/Rj  -  0.5  and  0.6  and  between  —  ■  0.2  and  0.4  ao  that  when  the 

eigenvalue  la  tracked  at  R^/Rj  *  0.5  froa  “  1  down  to  0.2  (figure  7)  and  then  at 

2 

—  -  0.2  froa  Rj/Rj  “  0.5  to  0.9  (figure  8)  we  arrive  at  a  different  eigenvalue  than  if 


we  follow  it  at  R^/Rj  "  0.9  froa  —  “  1  down  to  0.2.  Tracking  the  interfacial  eigen¬ 
value  for  Rj/R2  >  0.6  aa  in  figure  9  ylelda  valuee  in  column  c2  in  figure  8. 
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Figure  7 


*1/R2  “  °*5  '  *•  “  100  '  aR2  ”  1  '  0 

The  first  colan  track*  th*  lntarfaclal  eigenvalue  c  teem  Its  Talus  at  v  ,/U 2  ”  t  down 
to  U1/»i2  -  0.2. 

Th*  sscond  oolusui  tracks  a  ooa-f  laid  eigenvalue  cgG  racordsd  by  Sal  wan  and  Grosch  for 


Truncation  of  Chebyshev  polynomial*  Is  at  N 


19 


“l 

g2 

c 

C8G 

1.0 

.75  ,  0 

.81160, -.28413 

.98 

.75272, -.288821-3 

.81567, -.28105 

.96 

.75556. -.598088-3 

.81988, -.27791 

.94 

.75851, -.929718-3 

.82322, -.27471 

.92 

.76159,-. 128598-2 

.82872, -.27146 

.9 

.76481,-. 166928-2 

.83337, -.26815 

.88 

.76817, -.208228-2 

.83819,-. 264 78 

.86 

.77169, -.252838-2 

.84319, -.26136 

.84 

.77538,-. 301 108-2 

.84838,-. 25788 

.8 

.78330, -.410338-2 

.85938, -.250 78 

.6 

.83885,-. 141168-1 

.93110, -.21345 

.52 

.87350, -.2282 58-1 

.971 18,-. 19946 

.5 

.88404, -.259208-1 

.98254,-. 19639 

.48 

.89557, -.295678-1 

.99445,-. 19362 

.46 

.90825, -.339188-1 

1.0068  ,-.19123 

.44 

.92234, -.391868-1 

1.0196  ,-.18933 

.42 

.93815, -.456868-1 

1.0322  ,-.18807 

.4 

.95619, -.539128-1 

1.0440  ,-.18751 

.38 

.97740, -.647038-1 

1.0531  ,-.18740 

.36 

1.0040  ,-.796468-1 

1.0551  ,-.18568 

.34 

1.0456  ,-.994608-1 

1.0389  ,-.17766 

.32 

1.1013  ,-.985108-1 

1.0100  ,-.18494 

.3 

1.1483  ,-.909438-1 

.99625, -.19537 

.28 

1.1953  ,-.843918-1 

.99031, -.20426 

.26 

1.246  ,-.789048-1 

.98896, -.21315 

.25 

1.2734  ,-.765488-1 

.98945, -.21787 

.24 

1.3027  ,-.744628-1 

.99050, -.22279 

.22 

1.3682  ,-.711928-1 

.99400, -.23322 

.2 

1.4458  ,-.693738-1 

.99928, -.24428 

^Convergence  was  tested  at  P  ^  Ai  ^  "  0.34. 


Figure  8 


V 

U 


1 

2 


0.2  , 


Re  -  100 


,  oRj 


Node  0 


The  first  column  c,  above  the  behaviour  of  the  eigenvalue  which  was  tracked  from  the 
Interfacial  eigenvalue  in  Figure  7  for  R1/R2  “  0.5. 

For  Rf/Bj  *  °‘6,  th#  eigenvalue  which  is  tracked  froa  the  Interfacial  eigenvalue,  for 
example  in  Figure  7  and  9,  lie  in  column  Cj. 


*1^*2 

C1 

e2 

0.5 

1.4458, -.693738-1 

.99928, -.24428 

0.52 

1.4632, -.773848-1 

.98996, -.24079 

0.54 

1.4788, -.864978-1 

.98087, -.23635 

.6 

1.515  ,-.12375 

.95413, -.21492 

.65 

1.5391,-. 17317 

.92814,-. 18215 

.68 

1.5576, -.21326 

.90574,-. 15303 

.7 

1.5747, -.24298 

.88523,-. 13000 

.72 

1.5966, -.27278 

.85914,-. 10415 

.75 

1.6388, -.31302 

.80715,-. 641698-1 

.76 

1.6551, -.32463 

.78624, -.510108-1 

.78 

1.6906,— .34510 

.73897, -.257708-1 

.8 

1.7298, -.36252 

.68444, -.278178-2 

.84 

1.8206, -.39118 

.55456,  .326778-1 

.88 

1.9309, -.41243 

.40176,  .477008-1 

.9 

1.9933, -.41842 

.32101,  .456198-1 
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rigor*  9 


Jtf/Rj  -  0.7*  ,  H*  “  100  ,  alj  ■  1  ,  Nod*  0 


th*  first  colimm  track*  th*  lnt*rf*ei*l  aiganTaloa  e  fra*  it*  vala*  *t  u1  Ai  j  “  1  down 
t»  *  <*•*• 

Th*  second  colau  track*  a  ona-f  laid  alganralua  oM  r*cord*d  by  Calmn  and  Groach  for 
U1^2  ‘  ’• 
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*1 

**2 

c 

8 

o 

i 

i 

1.0 

.39160,  0 

.61160, -.26413 

.9 

.40351,  .143988-2 

.66443, -.26774 

.6 

.41611,  .323931-2 

.93120, -.23273 

.7 

.43647,  .336471-2 

1.0161  ,-.24078 

.6 

.46033,  .776331-2 

1.1266  ,-.23714 

.5 

.49279,  .101669-1 

1.2929  ,-.26247 

.4 

.33934,  .112711-1 

1.3296  ,-.33730 

.33 

.37139,  .101911-1 

1.3926  ,-.36369 

.3 

.61237,  .637601-2 

1.4673  ,-.36010 

.2 

.73897,  .237701-1 

1.6906  ,-.34510 

R1/R2  “  0.6  ,  *a  -  100  ,  oRj  -  1  ,  Nod*  0 

Tha  intar facial  alganralua  o  i*  trackad  from  it*  ralu*  at  v^/U2  m  1  down  to  0.2. 
Th*  ona-fluld  alganvalue  ou  i*  trackad  cloaaly  batman  "  0 .46  and  0.2. 


M1 

U2 

0 

e8G 

1.0 

.64  ,  0 

.8 

.67321, -.227106-2 

.6 

.73259, -.944791-2 

.92385, -.52808 

.46 

.7997 3, -.267496-1 

1.0639  ,-.40287 

.44 

.81229, -.314726-1 

1.0907  ,-.37493 

.42 

.82378,-. 372178-1 

1.1173  ,-.33474 

.4 

.84021, -.442678-1 

1.1142  ,-.29304 

.38 

.85338, -.529938-1 

1.1175  ,-.27473 

.35 

.87994, -.703098-1 

1.1377  ,-.23073 

.33 

.89628, -.854888-1 

1.1584  ,-.23392 

.3 

.91806,-. 11425 

1.2022  ,-.20574 

.28 

.92884,-. 13602 

1.2423  ,-.18576 

.23 

.93956,-. 16804 

1.3218  ,-.15788 

.23 

.94509,-. 18766 

1.3887  ,-.14216 

.2 

.95413, -.21492 

1.315  ,-.12375 
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*••1  (c)  dtcrMMi  with  a  and  approach**  »(*,)  aa 

(aaa  graph  12).  Only  for  R,/Rj  >  0.7  waa  tha  configuration  of  ainiaia  dlaalpation 
found  to  ba  atabla.  At  Ha  -  100.  a*2  -  1,  compare  R^Rj  -  0.7,  0.0,  0.9  (graph*  13,  6, 
14).  All  aodaa  hove  tha  aaa*  behaviour  for  R^/Rj  “  0.0,  0.9»  thay  ara  unatabla  If 
l>1/b2  <  1  and  atabla  If  V^Ata  >  1.  However,  w*  find  that  for  w  1/^*2  <  1  tha  aaroth  aoda 
can  baaaaa  atabla,  (aaa  Ha  ■  100,  R^/Rj  ”  0.2,  0.5  or  H1/R2  “  0.8,  grapha  15,  16,6). 

Largo  aR^,  moderate  Ha,  W1/b2  0  (0.2,  21 

Aa  aR2  lnoraaaaa,  tha  ragion  of  atabillty  dacraaaaa. 

Grapha  17  -  is  llluetrat*  tha  variation  with  R^/Rj  at  Ha  -  1,  aRj  -  10.  Grapha  19 
-  31,  illoatrata  tha  variation  with  R,/R2  at  Ho  -  100,  082  -  10.  In  practlca,  aurfaoa 
tanalon  any  daapan  tha  lnatabllitlaa  aa  in  tha  eaaa  of  flow  down  an  lncllnad  plana  (Tib,  p. 
502).  Por  aaall  a,  Rick ox  ahowod  that  aurfaoa  tanalon  affaota  ara  nagllglbla  oxcopt  whan 
R^/*2  <  0.1.  tha  raal  (c)  for  all  aodaa  lia  eloaa  togathar.  So*  Graph  22. 

Laroar  Ha.  At  Ha  -  1000,  R^/R^  “  0.8,  ORj  -  1,  graph  23  (coapar*  with  graph  6)  ahowa 
that  atabillty  la  loat  via  tha  0th  aoda. 

Larger  u 1/u2  tha  variation  of  la  (o)  with  larger  than  2  for  given  a  and 

Rj/Rj  la  laaa  rapid  than  for  aaallar  11  ^  Aij.  8oo  grapha  24-26  for  Re  “  100, 

oRj  -  1,  8,/Rj  "  0.5,  and  Ha  -  100,  oRj  -  1,  R^Rj  -  0.9  and  Ra  -  100,  oRj  -  10, 

Rj/Rj  -  0.9. 

tha  axparlaanta  reported  by  Rverage  (2J  are  for  M  2  *  6  and  equal  volume  flux**, 
fraa  which  wa  ccaputa  R^/Rj  -  0.57.  ihla  value  alao  agraoa  with  hla  photograph. 
Caloulatlona  in  Graph  24  ahow  that  for  Wj/Wj  >  4,  la  c  change*  very  little.  It  follow* 
froa  our  calculation*  (graph  24)  that  the  concentric  configuration,  in  the  experiment  of 
Kveraga,  ahould  ba  unatabla.  However,  the  growth  rata*  of  tha  unatabla  aodaa  ara  very 
aaall,  eo  lnatabllitlaa  aay  taka  a  longer  tiaa  to  develop  than  tha  observation  tla*  in 
Xveraga'a  experlaant.  It  la  alao  poaaibl*  that  the  inatability  only  grow*  to  a  aaall 
aaplltuda,  and  tha  eventual  flow  alght  not  b*  vary  different  from  tha  concentric  pattern. 


Zb  tht  mm  thin  la  larga>  **  would  ba  tha  mm  In  *  lubrication*  applications 

auoh  aa  oil  tranaport,  wa  find  that  concantric  aBMpaulatlM  la  tha  prafarrad 
configuration. 
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GRAPH  4 


0 


Imaginary  part  of  c  versus  viscosity  ratio  for 
low  Reynolds  number  and  moderate  a. 

Re  -  1,  aR^  -  1,  -  0.8 


Modes  0  to  4  increase  as  u./y,  decreases  from  1  to  0.2 


GRAPH  7 


Re  -  100,  aRj  -  1,  node  0 


Numerals  next  to  the  curves  denote  V^/u2 • 


GRAPH  8 


Re  ■  100,  aRj  "  1,  mode  1 


Numbers  next  to  curves  denote  — 

U2 


V*2 
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GRAPH  9 

Re  *  100,  aR2  -  1.  “°de  2 
Numerals  next  to  curves  denote  ^  . 

for  Hl.i.e  to  2.0  are  indistinguishable 
^2 

under  these  scales. 


GRAPH  12 

Real  part  of  c  versus  viscosity  ratio. 
10,  aR2  ■  1,  Rj/Rj  “  0.6,  velocity  scale 


GRAPH  13 


*•  ■  100,  aRj  -  1,  I^/Rj  -  0.7 

The  magnification  on  the  right 

W, 

displays  behaviour  for  1  <  —  < 


GRAPH  14 


Ra  ■  100,  Rj/Rj  “  0.9,  aRj  “  1,  velocity  se«le  *  1 
Curve  Cor  node  3  liee  between  modes  2  and  4 


GRAPH  22 

Real  part  of  c  versus  viscosity  ratio 
Re  “  100,  Rj/Rj  -  0.6,  oR2  «  10,  velocity  scale  «  1 

Curves  for  nodes  1,  2  and  3  lie  between  nodes  0  and  4 


0.8 


GRAPH  23 

Imaginary  part  of  c  versus  viscosity  ratio 
Re  «  1000,  aR2  ■  1,  Rj/R2  “  velocity  scale 
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GRAPH  26 


Ra  »  100,  R^/Rj  -  0.9,  aRj  -  10 
"Largo  a"  and  "large  y^/y^"  ragi 


Mod  a  a  0  to  4  hava  tha  iut  bahavlour 
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